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o 
o 

o 

Q 



Maltc Vogl, Gemot Schallcr,* and Tobias Brandes 
Institut fur Theoretische Physik, Technische Universitdt Berlin, Hardenbergstr. 36, 10623 Berlin, Germany 

For sufficiently low reservoir temperatures, it is known that open quantum systems subject to 
decoherent interactions with the reservoir relax towards their ground state in the weak coupling limit. 
Within the framework of quantum master equations, this is formalized by the Born-Markov-secular 
(BMS) approximation, where one obtains the system Gibbs state with the reservoir temperature as 
a stationary state. When the solution to some problem is encoded in the (isolated) ground state of 
a system Hamiltonian, decoherence can therefore be exploited for computation. The computational 
complexity is then given by the scaling of the relaxation time with the system size n. 

We study the relajcation behavior for local and non-local Hamiltonians that are coupled dissipa- 
tively with local and non-local operators to a bosonic bath in thermal equilibrium. We find that 
relaxation is generally more efficient when coherences of the density matrix in the system energy 
eigenbasis are taken into account. In addition, the relaxation speed strongly depends on the matrix 
elements of the coupling operators between initial state and ground state. 

We show that Dicke superradiance is a special case of our relaxation models and can thus be 
understood as a coherence-assisted relaxation speedup. 

PACS numbers: 03.65.Yz, 31.15.xp, 03.67.-a 
Keywords: Lindblad form, stationary state, dissipation 
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I. INTRODUCTION 



Since the seminal papers by Shor [1] and Grover [2] 
the interest in quantum computation is unbroken and 
has given rise to plenty of research [3] . 

However, the superiority of quantum computation re- 
lies on the fragile quantum coherence of states in the 
system of interest, which is sensible to interactions with 
the environment [4]. The process of loosing coherence, 
the so called decoherence, has yet to be fully under- 
stood [5, 6]. Various attempts are constantly made to 
describe the process of decoherence for special systems 
analytically [7]. The hope is that one can overcome the 
devastating influence of decoherence by error correction 
codes [8] - or by developing new models of computa- 
tion such as measurement-based [9, 10], holonomic [11], 
ground-state [12, 13] or adiabatic [14] quantum compu- 
tation. 

In this paper we describe the (unitary and non-unitary) 
effects of decoherence by means of a Markovian master 
equation, namely the Born-Markov-secular (BMS) ap- 
proach [15]. The resulting master equation is an effec- 
tive evolution equation for the reduced density matrix 
of the system. Such master equations arc valid in the 
weak-coupling and long-time limit provided the reservoir 
relaxation time is significantly faster than the system dy- 
namics. The BMS master equation is always of Lindblad 
type [16] and thus preserves positivity of the reduced 
system density matrix at all times. Beyond this it can 
be shown, that for a thermalized reservoir with reservoir 



temperature /3 the equilibrated system Gibbs state 



-PHs 



Pth 
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(1) 



with the same temperature /3 is a stationary state of the 
BMS master equation [15]. When we label the energies 
of the system Hamiltonian with a unique ground state 
hy Eq < El < . . . < Em-1 for an A^-dimensional system 
Hilbert space, the probability to find the stationary state 
in the systems ground state is given by 



Pn 



> 



1 
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where the lower bound is a rather crude estimate. As- 
suming a system composed of n qubits (such that N = 
2"), this implies a large probability to reach the ground 
state by dissipative relaxation when the temperature 
/cbT ~ (3^-^ scales as 



k^T{n) oc 



[El - Eo) 



(3) 
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This remarkably mild dependence of the required tem- 
perature on the system size n is not spoilt by the fun- 
damental energy gap Ei — Eq: In fact, it is possible to 
encode the solution to interesting problems in the ground 
state of local Hamiltonians with a fundamental energy 
gap that does not even scale with the system size n at all. 
Since such Hamiltonians can easily be constructed for the 
NP-completc problems Exact Cover 3 [14, 17] and 3-SAT 
[18], it does not come as a surprise that also the solution 
to the prominent factoring problem can be encoded in 
the ground state of a two-body local Hamiltonian [19]. 

The robustness of the ground state against decoherence 
at low temperatures is already exploited in various com- 
putation schemes (see e.g. [12, 13, 14, 20, 21]). However, 
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here we are aiming at solving problems using nothing 
but decoherence-induced relaxation [22]. Given a corre- 
sponding system Hamiltonian, the question for the speed 
of this relaxation process arises naturally. The depen- 
dence of the relaxation speed on the system size defines 
the computational complexity of this cooling algorithm. 

In this paper we investigate the scaling behaviour of 
the relaxation speed for several time-independent (and 
scalable) Hamiltonians. The paper is organized as fol- 
lows: In section II we introduce our methods, which we 
apply to nonlocal and local Hamiltonians in sections III 
and IV; respectively. 



II. METHOD 

We separate the Hamiltonians for the various problems 
in this paper as H ~ Hg + i?sB + j where Hg and 
only act on the system and bath Hilbert space, respec- 
tively. In contrast, the interaction Hamiltonian can be 
decomposed as 



(4) 



A 



with the small dimensionless coupling parameter A ^ 1 
and the system {A_a) and bath (BX) coupling operators. 
By a suitable redefinition [15] we can always assume her- 



mitian coupling operators Aj^ = and Bj^^ 
define the bath-correlation functions as 



b\. We 



CABir) = TrB {e*^-^i?^e-*^--BepS^} = , (5) 

and introduce their even and odd Fourier transforms 
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where V denotes the Cauchy principal value. For the 
transformation into the interaction picture, it is useful 
to expand the system coupling operators into eigcnopcr- 
ators of the system Hamiltonian 



A. 



AAiuj) 



A 



(7) 



which obey by definition [ffg, ^^(w)] = — where 
ui denotes all energy differences of the system Hamilto- 
nian. Such a decomposition can for example be achieved 
in the system energy eigenbasis Hs \a) ~ Ea\a) by using 



(a I Aa \ b) {b\. Performing 



the Born-Markov and the secular approximation, where 
the latter amounts to the neglect of the resulting oscil- 
lating terms [15], we arrive at a Lindblad form master 
equation for the system density matrix 
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where cr^B(a;) is anti-hermitian (leading to a hermitian 
Lamb-shift Hamiltonian) and 7^6(0;) is positive semidef- 
initc for each to. The Lindblad form of this equation is 
a sufficient condition [15, 16] for preserving positivity of 
the density matrix for all times and coupling constants 
A. When we phrase these results in the system energy 
eigenbasis, the quantum master equation takes the form 
(with pab={a\p\b)) [23] 
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where the damping coefficients equate to 

7ab,cd = ^"^^lABiEb ^ Ea)5Ed-E^,Ei-E^ 
AB 

X {a\AB\b) {c\AA\d)* 
and the Lamb-shift terms are explicitly given as 



,(9) 
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We will refer to Eq. (9) as the quantum master equation 
in the following. 



Dephasing model 

As a simple example for the quantum master equation, 
let us consider a pure dephasing model, where the cou- 
pling operator commutes with the system Hamiltonian. 
This implies that Hg and A can be diagonalized simul- 
taneously A\a) = da I a), and the coefficients acquire an 
especially simple form 



lab,cd = >^^l{^)5ab5cddadc , 



&ab = ^^(0)<5a6da ,(12) 
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such that we obtain for the matrix elements of the density 
matrix 

I \{Ea - Eb) + ^a(O) (dl - 4) 



Pab 



2i 



A^(0) 



[da 



Pab ■ 



(13) 



which leads to an exponential decay of off-diagonal ma- 
trix elements that correspond to different eigenvalues of 
the coupling operator A (recall that cr(0) is imaginary) 



Pab{t) oc exp 
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Rate Equation 
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operators to avoid non-ergodic behavior (many station- 
ary states) [15]. Regarding the system coupling opera- 
tors, we note that as p{t) — e'~* p'^ and C cx a scaling 
of the coupling operators has a direct effect on the time- 
dependence of the density matrix. In order to perform a 
fair comparison between models with different coupling 
operators, we therefore concentrate - unless noted oth- 
erwise - on operators with a strictly bounded spectrum 
(i.e. ||j4_4|| < 1, where the spectral norm for a hermitian 
operator is defined as ||A^|| = v'^max)- the following, 
we will consider models with a unique ground state de- 
noted by \w) throughout the paper. Furthermore, to be 
able to compare numerical with analytical solutions, we 
define the relaxation time trciax for a given initial state 
as the time needed, to achieve a ground state population 
of Pground = {w\ p\'w) = 0.9 in thc stationary state p. 



It is evident from Eq. (9) that for a non-degenerate 
spectrum of the system Hamiltonian the BMS dynamics 
of the diagonals (populations) 
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completely decouples from the dynamics of the off- 
diagonals (coherences) and is also unaffected by the 
Lamb-shift terms. Therefore - compatible initial con- 
ditions provided ~ the coherences may be completely ne- 
glected from the considerations. 

In reality, it can be expected that exact degeneracies 
of ideal model Hamiltonians may be broken by pertur- 
bations. In addition, as shown before, dephasing noise 
may also lead to the decay of coherences, such that the 
above equation can be well motivated. We will refer to 
Eq. (15) as rate equation in the following. Note however, 
that even for a non-degenerate spectrum with smallest 
level splitting Ae, the secular approximation is only ap- 
plicable for times t > Ae~^, which should be kept as 
an additional overhead in mind when discussing the scal- 
ing of relaxation. For a non-degenerate two-level system, 
rate equation and quantum master equation coincide. 



General Assumptions 

Throughout the paper, we will assume a bosonic bath 
Hb — J2k (b\bk + constantly held at thermal equi- 
librium. For simplicity, we will consider single-operator 



hkbk + Kbl 



such that we 



couphng = XA®J2k 
may omit the indices and the Fourier transform of the 
bath correlation function is given by [23] 



7-46 (^) = = 



.9(1^1) 



(16) 



[l-e-'3"[ ' 

where we have assumed a thermalized bath with inverse 
temperature /?. It should be noted that for large sys- 
tems, one will in general need many different coupling 



III. NONLOCAL COUPLING OPERATORS 

All coupling operators in this section will - an imple- 
mentation based on n qubits assumed - be nonlocal in 
the sense that they act on all qubits simultaneously. In 
addition, most of the following problems will have the 
symmetry that the dynamics is completely contained in 
the subspace spanned by [ w) and [ s) , where 



1 



)(|o) + |i)),^:^E 



N 



N 



(17) 



is the superposition of all computational basis states, and 
\w) is the distinguished ground state. An orthonormal 
basis for this subspace is given by 



\V2) = 



(wis) Iw) 



= Ei«)-(i8) 

a^w 



We will consider a system Hamiltonian with the unique 
ground state \w) and exponentially many degenerate ex- 
cited states, which in basis (18) can be written as 



Hs ^AE[1- \w) {w\] 




AE 



(19) 



where \w) may for example be a marked item in a 
database [24] or encode the solution to a problem, such 
that the Hamiltonian may act as an oracle [2]. A classi- 
cal minimization algorithm would find the solution \w) in 
the worst case after N trials of this oracle. Accordingly, 
we have £"0 = and Ei = . . . = E^_i = AE > 0. In 
order to get a large overlap between the Gibbs state (1) 
and the ground state, we require low temperatures 



BAE 



OiN)., 



(20) 



where N = 2" denotes the dimension of the Hilbert 
space, such that the required temperature must scale in- 
versely with the system size, compare Eqn. (3). 
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A. A single Projector 

We consider a single coupling operator with an inter- 
action Hamiltonian of the form 



A 



1 / 1 ^/N - 1 



TV V \/N -1 N -1 



(21) 



As it is a projector, the operator A has eigenvalues zero 
and one. We introduce the rate equation variables 



Pu 



(22) 



and the quantum master equation variables variables 

(23) 



qc 

^1 = Pww 1 



'2 = 5Z P°-^ ■ 
a, 6^10 



with which we can both write Eqs. (15) and (9) in the 
form 



Z2 



Zl 



(24) 



where M is a two by two matrix. Note that trace con- 
servation implies that z^-^ + ^2° = Oi but -I- is not 
conserved. The formal solution of Eq. (24) is given by 



(25) 



where denotes the initial conditions. As the solution is 
unknown, only one initial condition that does not favor 
the ground state can be prepared in the rate equation 
approach. Here 



1 

TV ' 



N ~ 1 
TV 



(26) 



which corresponds to a diagonal density matrix with 
equal diagonal elements - or simply the Gibbs state (1) at 
infinite temperatures. In the quantum master equation 
the additional initial state 



1 

TV ' 



TV 



(27) 



is also preparable. This state corresponds to a |s) (s| 
initial density matrix, as can be seen from Eq. (23). 



1. Rate Equations 

For the rate equations we obtain with variables (22) 
from Eq. (15) the matrix 



A^/-7(-Ai?)(iV-l) liAE)\ 

7V2 I 7(-Ai?)(TV-l) ~j{AE) ' ^^^^ 



Without knowledge of the solution \w), we can prepare 
the initial condition (26). Solving Eq. (25) with the ma- 
trix (28) we find 



zi{t) 



1 



(TV + e 



x(e'3^^-l)e- 



I3AE 
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1) TV 

, (JV + e/3^E-l) 



[(TV-1) 



(29) 



such that the steady state (1) will be reached after time 
t 3> Tio, where 



TV2 fe/3A£ _ l^ 



X^g{AE) (TV-t-e/sAB- 1) 



O (TV^) , (30) 



as can be seen from the exponent. This is surprisingly in- 
efficient. In addition, it must be kept in mind that when 
we assume the exact degeneracies of Eq. (21) to be lifted 
by some perturbation, the rate equations become valid 
after times larger than the largest inverse level splitting, 
which would lead to an additional computational over- 
head. This overhead will be significant, as the splitting 
must be exponentially small: The 2" — 1 non-degenerate 
excited states must fit into the range [AE — 5, AE + S] 
with S <^ AE when the Hamiltonian (19) should be ap- 
proximately valid. However, alternatively we may as- 
sume that dephasing noise leads to the decay of coher- 
ences and thereby motivate the simplified rate equations 
after a significantly shorter timescale. 



2. Quantum Master Equations 

With the coherences, relaxation is more efficient. With 
variables (23) we obtain from Eq. (9) the evolution ma- 
trix 



(M ^^2 / -7(-Ag) +2(±A£)\ 

- 1U2 +7(-A-E) -7(+A-E) ' l-^-^i 

\ (1) (N-l) ) 



see also Appendix A. For both initial conditions (26) 
and (27), we obtain a decay into the steady state for 
times 



iV2tanh(^) 
\^g{AE){N-\) ^ 



O (N) , (32) 



which is not faster than classical minimization. How- 
ever, for the previous initial condition (26), the overlap 
between steady state and ground state becomes expo- 
nentially small. This can be understood since the first 
choice (26) is not fully contained in the subspace spanned 
by the dyadic product of (18), such that only an exponen- 
tially small part of this initial condition will relax towards 
the Gibbs state. In contrast, the initial condition (27) 
corresponds to \s) {s\ and is thus fully contained within 
the relevant subspace. In this case, relaxation approaches 
the Gibbs state (1). This problem of the evolution of 
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the quantum master equation (9) towards a non-thermal 
state can be cured by using couphng operators that medi- 
ate transitions out of the two-dimensional subspace [15]. 
Note that the evolution matrices for both methods (28) 
and (31) coincide for N = 2, that is for one qubit n = I. 



D. Hadamard Coupling 

Inspired by the nonlocal representation of the previous 
coupling operator (21) we choose a product of Hadamard 
gates as coupling operator 



A 



B. An indirect channel 

As a slight modification of the previous coupling oper- 
ator, we consider 



A = Tj{\s) {s\ + \w) (w\) 



(33) 



In order to have ||A|| < 1, we may choose r] — jj^, which 
converges to one for large N. 

With variables (22) Eq. (15) leads to the evolution ma- 
trix M = ■jjJ^^M"', where M''° is given in (28). There- 
fore, with initial conditions (26), the rate equation relax- 
ation time will scale as t ^ 0{N'^), which is the same as 
in Eq. (28). 

With variables (23) the evolution equation correspond- 
ing to the quantum master equation (9) can be written 
as M = ■p^^^Af™'', with A/'°° as given in Eq. (31). 
Again, we have the same scaling of the relaxation time 
as in Eq. (31) - 0{N). The dependences of the relaxation 
time on the number of states for both cases are depicted 
in Fig. 1. 



C. A direct channel 

Here we will demonstrate that the runtime can become 
independent on the system size when the coupling opera- 
tor mediates a direct transition between the initial state 
and the ground state. In this section we focus on this by 
considering the coupling operator 



V2 



mi 



(35) 



The eigenvalues of the operator A are given by —1 and 
4-1 (with exponentially large degeneracies) and its matrix 
elements in the computational basis read 



{a\A\b) = n 



(-l)fci _ (_l)ao6 



N 



N 



(36) 



where a o b = X]a=i '^aba denotes the Hamming weight 
(the number of ones in the binary decomposition) of the 
bitwise product of a and b. 



1. Rate Equations 

In analogy to Subsec. Ill, we obtain due to the different 
matrix elements |(a| A |6) | ^ 1/N with the variables (22) 
the evolution matrix M = NM"", cf. Eq. (28). For the 
initial equal distribution (26) we obtain relaxation into 
the steady state (1) after 



t > 



N 



OiN) , 



(37) 



with Tro given in Eq. (30), which corresponds to the classi- 
cal complexity. Note that the steady state is thus reached 
in similar time as for the quantum case of Subsec. Ill A, 
as can be seen from Fig. 1. 



2. Quantum Master Equations 



A = Tj{\s) (w\ + \w) (s\) 



(34) 



where it is obvious that (.s| A jw) = Ojl}, since in order 
to have ||A|| < 1, one may choose the prefactor rj = 

1+Vn' 

Using the variables (22) we can write Eq. (15) with 

rM'°. Therefore, the 



N 



the evolution matrix AI — — t= — tttj 

(VA'-l-l)'' 

scaling of the rate equation dissipation time is given by 

With the variables (23) we obtain for Eq. (9) the evo- 
lution matrix M = —^^——Ad™". For initialization with 

an \s) {s\ density matrix as in Eq. (27) one obtains a large 
overlap of the ground state and the stationary state, and 
the relaxation time scales as t ~ 0(1)! These scaling 
behaviors are summarized in Fig. 1. 



Due to the different matrix elements of our coupling, 
we cannot close Eq. (9) with the variables (23). Instead, 
one achieves closure with the variables 



zi = Pu 



Z2 



Pab- 



(38) 



ab^w 



Insertion into Eq. (9) yields an evolution equation of type 
z = AIz with AI = AfM"^<=, cf. Eq. (31) and see Ap- 
pendix B for details. Again, we do now have different 
possibilities of initialization without knowing the actual 
ground state w. 

For example, by using a diagonal density matrix as 
initial condition as in Eq. (26) we obtain relaxation after 
times 



t»^^'^^0(l) 



(39) 
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with Tme as given in Eq. (32), but the stationary state 
always has exponentially small overlap with the ground 
state. 

Even when we initialize the density matrix with |s) (s| 
(see Appendix C for the details), the overlap of the sta- 
tionary state with the ground state \w) will be expo- 
nentially small - unless the ground state is given by 
\w) = |0) (g) . . . (8) |0). In the latter case, relaxation 
also becomes extremely fast, since the relaxation time 
becomes independent of the system size. This is due to 
the fact that the Hadamard coupling provides a direct 
channel from the initial state towards the ground state, 
as \s) =H^"|0...0). 
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FIG. 1: [Color Online] Scaling of the runtime to achieve 
Aground ~ 0.9 versus the number of states A'^. For each A^, the 
inverse temperature f3 has been adapted such that the ground 
state probability of the Gibbs state (1) was Pground ~ 0.95. 
For all couplings one observes a clear speedup for the master 
equation results (hollow symbols) in comparison to the rate 
equation results (symbols) . The rate equation solution for the 
Hadamard coupling (blue ^-symbols) is of the same order as 
the quantum master equation solution of Subsec. Ill A 2. The 
runtime for the quantum master case the Hadamard coupling 
(blue triangles up) is only finite if |0---0) is chosen as the 
solution. By construction, the quantum master equation run- 
time is always finite for the operator of Sec. HI C (green dia- 
monds). Dashed lines denote the asymptotic scalings T cc N 
and T cc N"^ , respectively. Other parameters have been cho- 
sen as fl(|A£;|) = 2, A£ = 1, and A = 0.01. 



IV. LOCAL COUPLING OPERATORS 



where Tia denotes the subspace of all states with Ham- 
ming distance a to the solution \w). The Hamming dis- 
tance is a relative distance measure and is defined as 
the number of bitflips needed to transform one com- 
putational basis state to another. Obviously, it ranges 
between (equal states) and n (inverted states). The 
above coupling operator mediates transitions between 
states with Hamming distance one. In order to constrain 
the eigenvalues of A between —1 and 4-1, we need to 
choose 77 = i. For the ease of manipulation, we or- 
der all states in the computational basis with respect 
to their Hamming distance to the solution, such that we 
can assume that Ei^ < Ei < . . . < En-i < En- There are 

states with Hamming distance a to the solution 

which also gives the degeneracy of Ea- In the fol- 
lowing, we will use that for each state in Ti^, there are 
(n — a) single-bitflips that lead to Ti-a+i and [a) single- 
bitflips that lead to Ha-i- 

As an example of this, consider n = 10 spins and the 
ground state 1 0000000000). Then, the spin configuration 
|1011010101)e7^6- There are six bitflips leading to 7^5, 
namely all bitflips on ones. In addition, the four possible 
bitflips on zeros lead to 7^7 . This generalizes to arbitrary 
solution states, as the Hamming distance is a relative 
measure. 

Note that already the simplest classical algorithm for 
finding the ground state is quite effective, as one only 
has to flip single bits of an n-bit bitstring and calculate 
the corresponding energy: If the energy decreases, one 
keeps the modified bitstring configuration, whereas in the 
opposite case one tries to flip another bit. Thus, one finds 
the ground state classically after at most n bit flips. 

Assuming an equidistant level spacing Ea ~ aAE, we 
see that the condition on the temperature to yield a large 
overlap between the ground state and the Gibbs state is 
significantly weaker 



0{n} 



(41) 



than in Eq. (20). 



A. Climbing a Hamming ladder 

Here we consider a Hamiltonian that gives an energy 
penalty Ea to all states that have Hamming distance a 
to the solution \uj) 



1. Rate Equations 



We introduce the variables 



= E E i«) («i ' ^ = ' (40) 



e=i 



E 



(42) 
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Then, we obtain from Eq. (15) 



time, the time necessary to reach a population of 0.9 in 
the ground state, for up to 400 qubits. The comparison 
as depicted in Fig. 2 shows the clear speed-up from 0{n'^) 
to 0{n) for the quantum case. 



+X^7^\n -a + l)7(£;„„i - E^)z'^'_ 
-X^if [{a)^{E^~E^.i) 
+ (n-a)7(i;„-£;„+i)]z™ 

+X^if{l)^{E,,^,-E,,)z^:_^ 
-XW[{nh{En-E,,.,)]zi:. 



(43) 



When we consider equidistant spacings — i = AE, 
the eigenvalues of the associated matrix are given by 



A. 



aX^ifg{AE) coth 



(44) 



where a S {0, 1, 2, . . . , n} such that we obtain an efficient 
(polynomial) scaling of the relaxation time t cx 7y~^ = th? 
in this case. 



2. Quantum Master Equations 
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FIG. 2: [Color Online] Scaling of the runtime to achieve 
^success = 0.9 versus the number of qubits n. For each n, 
the inverse temperature p has been adapted such that the 
ground state occupancy was Pground = 0.95. One observes a 
clear speedup for the quantum case (black squares) in com- 
parison to the rate equation case (red circles) , where only the 
populations are considered. Dashed lines denote the asymp- 
totic scalings T <x n and T oc , respectively. Other parame- 
ters have been chosen as (;(|Ai?|) = 2, AE — 1, and A = 0.01 



We introduce the variables 



(45) 



Then, we obtain from Eq. (9) 

zl'' = -^X^v^n-f{E^ - E,)zl'' + XWliEi - E^)zT 



+X'if{n - a + lYi{E^-i - E^)zl'Li 



-xW 



-/{Ea - Ea-i)a{n - a + \) 



+jiE^ - Ea,+i)in -a){a + 1) 
+XS^ia + ifjiE^+i - ^a)zr+i 

+A27y2^(S„_i - E,-,)zZi 
~X^jf^{E,,-En-i)zl^ . 



(46) 



To make a statement about the efficiency of this algo- 
rithm, we use numerical solutions of the coupled system 
of differential equations. However, due to the tridiago- 
nal structure of the system, standard solvers are quite 
efficient. Therefore, we can easily compute the success 



B. Example: Dicke Superradiance 

The question of relaxation speed has also been raised in 
the field of quantum optics, since it is connected with the 
energy flux radiated by the system into the reservoir. We 
have observed that with coherences the relaxation speed 
is quadratically enhanced. In the field of quantum op- 
tics, this is already known as Dicke superradiance, which 
exemplifies the broader validity of our approach. There, 
coupled relaxation equations of tridiagonal type also arise 
naturally: 

The Hamiltonian of Sec. IV is a generalization of the 
multi-mode Dicke Hamiltonian 

Wo 



i^Dickc = -j-r + x,p ®Y.^hkhl + hibk) 

k 



(47) 



where we have J" = ^f- The phenomenon of su- 

perradiance as first pointed out by Dicke [25] describes 
the process of coherent radiance by a system of n two- 
level systems, that are initially collectively excited. The 
radiance, contrary to classical emission, has an inten- 



sity / 



and a delay time tu before it reaches its 



8 



maximum [26]. The radiated energy current can be de- 
rived from the time-derivative of the expectation value 
for the energy. To observe these effects, wc initiahze 
the system in the state of highest Hamming-distance 
{z'^ = {0, • • • ,0, 1}). The original Dicke effect refers to 
zero temperature, i.e., we consider the limit /3 — > oo here. 
In this limit, we can easily solve Eqs. (43) and (46) ana- 
lytically by e.g. Laplace transform [27]. The solution of 
a system of the form 



yo 



PoVo + loVi , 



Vi = Pzyi + -fiVi- 



+1 ■ 



Pnyn 



(48) 



subject to the initial condition y° = (0, . . . , 0, 1) is then 
given by 



Vk 



(0 = E 



a—k 



n 

b—k:b^a 



Pa-f3b 



'n-1 
.c— 



,(49) 



where we have assumed nondegenerate coefficients Pa 7^ 
Pb^a- In case of Pa — Pb one may cither use higher-order 
residue formulae or simply analytic continuation of the 
above result. 



1. Rate Equations 

In the energy eigenbasis, we can easily relate the expec- 
tation value of the energy with the variables introduced 
in Eq. (42) via 



a=0 c^Hc 



(50) 



'ri Jh 

Ei^a4^(t)=t^oE("-f)<'^w- 



a=0 



From Eq. (43) wc identify in Eq. (49) in the 
zero-temperature limit Pa = —aX^r]^ g{ujQ) and 
7„ = +{a + l)\^7fg{u;a). 



2. Quantum Master Equations 

The quantum master equation expectation value of 
the energy can be obtained by noting that the Dicke 
Hamiltonian (47) is symmetric under all permutations 
of the spins, which implies that the whole evolution can 
be expressed with the angular momentum eigenstates 

|j, to) = 2to Ij, m), where j = n/2. These eigenstates 



relate to the superposition of all states in a given Ham- 
ming subspace via the relation 



n n 
2' 2 



1 




E 



(51) 



where a e {0, 1, 2, . . . , n} denotes the Hamming distance 
to the ground state |1, . . . , 1). For example, we see that 
for 71 = 4 the states 



|2,-2 
12 



[1111 

1 

71 



|2,+2) 



1 

10000) 



1110) 4- |1101) + |1011) + 10111)] 
0011) + 10101) + |1001) 
1100) + |1010) + (0110) , 

0001) + [0010) + [0100) + |1000)] 



(52) 



are symmetric under all permutations of the qubits. This 
can be exploited to relate the master equation expecta- 
tion value for the energy with the variables introduced in 
Eq. (45) via 



2 

n 



Tr{J>} 



Q=0 



E 



^oE 



a=0 



n\ In 
p\-,a 



2 I 

E 

a,beH^ 



2 

;a|p|6) 



'it) 



(53) 



From Eq. (46) we identify in Eq. (49) for zero temper- 
ature the coefficients Pa = —a{n — a + l)X'^ri^g{LUo) and 
7q = +{a + lY\^rfg{uLii)). These can be used to obtain 
the short-time dynamics from Eq. (49): Obviously, we 
n 
a 

for large a will contribute 



have < 



such that initially only the variables 



ZT + 



1 



(54) 



Performing the time derivative and assuming large n, we 
get an approximate expression for the radiation intensity 



-2^'9("°)"*g(a.o)A2n2c^o (e+^'f^"")"* - l) (55) 



For this approximate expression, we can analytically ex- 
tract the position and height of the maximum as well as 
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the width of the radiation peak as 



^peak 
-^pcak 



ln(2) 



cx 



— A g[ojQ)n (X n , 
2 In 



A 



peak 



V2+1 



1 

oc — 
n 



(56) 



If one considers the uppermost three variables instead of 
the uppermost two. one obtains similar scaling relations. 
These rough scaling relations for the Dicke flash are well 
reflected in the full solution, see Fig. 3, where also the 
solutions for the rate equation radiation intensity are dis- 
played for comparison. 




n = 20 
n = 40 
n = 80 



10' 

time t [a.u.] 



10 



10" 



FIG. 3: [Color Online] Plot of quantum master equation 
(solid) and rate equation (dashed) radiation intensity versus 
time for different system sizes n, with n = 20 (black thick 
line), n = 40 (red medium line) and n — 80 (blue thin line). 
Whereas in the incoherent case one always obtains simple 
exponential decay (dashed lines), the coherent case exhibits 
a strong flash with intensity proportional to and width 
proportional to 1/n after a delay time. The totally radiated 
energy (area below the curves) is the same for both coherent 
and incoherent emission. Other parameters have been chosen 
as f3ujo = 100.0, t^o = 1.0, A = 0.1, g{ujo) = 1-0. 



evolution was always ergodic. In order to restore er- 
godicity, more different coupling operators (decoherence 
channels) are generally needed [15]. Note that - slight 
energy splittings provided - ergodicity is automatically 
provided, since for times larger than the inverse level 
splittings the rate equation (with ergodic evolution) be- 
comes valid. A smooth transition between quantum mas- 
ter and rate equations is given by the dynamical coarse- 
graining method [23] or also found in the singular cou- 
pling hmit [28]. 

For the non-local coupling operators considered, we 
did not obtain the algorithmic efficiency of the Grover 
search. Note however, that this is possible in dissipa- 
tive approaches, either phenomenologically (i.e., without 
deriving the Liouville superoperator from the interaction 
Hamiltonian) [29] or also for an evolution under a Hamil- 
tonian [21]. In addition, not even the classical scaling 
was met by the rate equation approach for local opera- 
tors. However, this originates in our strict constraints on 
the coupling operator norm. 

Furthermore, we note that the algorithms can be ex- 
tremely efficient, when the matrix element of the coupling 
operator between the initial state and the ground state 
(s| is large. When these matrix elements become 

independent of the system size n, we may also observe 
relaxation times that become independent of n. These 
findings could be interesting for ground state quantum 
computation, where we can encode the solutions to diffi- 
cult problems in the ground state of a Hamiltonian. 

When considering local operators, we have re- 
discovered the Dicke superradiance Hamiltonian as a spe- 
cial case, such that one may now generalize the superra- 
diance not only to finite temperatures, but also towards 
non-equally distributed energy levels. Note however, that 
in the latter case the permutational symmetry of the 
Dicke Hamiltonian is destroyed, such that the radiance 
will then have to be derived by means of an n-resolved 
master equation approach [30]. 



APPENDIX A: SINGLE PROJECTOR: 
CALCULATION OF z FOR THE QUANTUM 
CASE 



V. CONCLUSION 

We observe a clearly faster relaxation speed when co- 
herences are considered, i.e., Markovian relaxation may 
be significantly slowed down by dephasing noise or - if 
degeneracies are slightly lifted - in the later stages. The 
heuristic explanation is, that coherences open more chan- 
nels of transitions for the system to reach the ground- 
state. In the simple examples we have considered, this 
speed up goes along with a non-ergodic evolution of the 
quantum master equation, i.e., the stationary state de- 
pends on the initial state, whereas for the rate equations 



From Eq. (9) we obtain 



Pww / ^cw^cwPww ~r y ^wc.wdPcd 

(Al) 

= -^(N - l)j{-AE)p.^^ + ^li+AE) J2 Pab , 

ah^w 



which is coupled to X^afc^tu Pab- Thus we also need the 
time-evolution for zi = J^at^w P"-b which can be calcu- 
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lated to 



Pab 



'Jca.db 2 '~ica,cd ''Icd.cb 



E E 

ab^w cd^w 

^ ^ ^aw , bio Pww 
ab^w 

{N-IMAE) ^ Pat 

ab^w 



Pab 



(A2) 



The Lamb shift terms cancel and the system of equations 
is closed, such that wc obtain Eq. (31). 



APPENDIX B: HADAMARD COUPLING: 
CALCULATION OF z FOR THE QUANTUM 
CASE 



Note that we make extensively use of 

1 ^(_^)aoc+6oc^^^^^ 



TV 
1 

TV 



V E^ ^ 



aoc+boc 



'5a(,-^(-l)"''°"', (Bl) 



whereas the time-evolution for Zi{t) is given by the fol- 
lowing, where aij and jij^ki were inserted: 



2i N LV ^ 



Pcb 



abc^w 
-(-l)'"°"+"'>ac] 

.A2a(0) 



2i N 

abcd^w 
^ -^doc-^-dob-ywoa-^-wob 



E\ { 1 \doa-\-doc^woa^xoob 
L(-i) Pcb 

PaA 



N 



J—L—L ^ ^ ^ ■^'■jWoa-\-wO: 

abcd^w 

1 2 7(+Aii;)(TV-l) 



b-\-aoc-\-bod 



Pcd 



-2' 



N 



Y^i-ir — Pcb 



2^ — (-1) 

abcd^w 



woc+wob 

bc^w 

-\-wob-\-cod+coa 



Pdb 



1.27(0) 
2 N 

l^^ -l{ + AE){N-l) ^ ,..^oa+»oc. 



ac^w 



y^^i Pad-yob) 

abcd^w 



which is derived from the property 



A 



2 A(^f + ^l)' 



(B2) 



This can be written as 



\-cob 



Therefore, one obtains 

J2{a\A\c){c\A\b) =^1(-1)-+^ 

C C 

= {a\A\b)=5ab- (B3) 
The time-evolution of zq ~ p^w is with Eq. (9) given by 



^ ^ ^cw,cw 



Pww ~i" ^ ^ ^wc,wdPcd 
cd^w 



-\^^{~AE) 



N ~ 1 



^^2 7(±^ ^(_l)™oa+.o.^^^^ 



ab^w 



Zl = 



-U^) E 



abc^w 



^ l^^iuoa+ujoc _|_ _|_ ^ -^ywob+woa 



Pab 



-X^-f{-AE) 



(TV -1)2 „7(+Ai?)(TV-l), 



+A 



2 7(0) 



TV 



TV 
1 1 

2^2 



-Zq - X 



Zl . 



N 



-Zl 



(B6) 



In the above equations, the Lamb-shift terms cancel and 
the system closes, such that we obtain the results of Sub- 
scc. HID 2. 
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APPENDIX C: NON-ERGODICITY: OVERLAP where d„ denotes the Hamming distance of the solution 
OF SOLUTION WITH THE STEADY STATE \cj) to the state |0 • • • 0). 

For an initial |s) (s| density matrix, the initial variable 
Z2 in Eq. (38) equates to 



aow~\-bow 

N ' ' 



a^iu l—l h^w k—1 

/ n 

n[ ■■■ n(-ir°"" - (-1)'" 

\aie{0,l} a„e{0,l}'=l 

U{N-lf for H = |0...0) 

(7V-2)4„A..o + ^, (CI) 
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